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Abstract 

f— ^ ■ The ^Li total photoabsorption cross section is calculated microscopically. As nucleon-nucleon 

00 ' 

^^ . interaction the semi-realistic central AV4' potential with S- and P-wave forces is taken. The inter- 

^. 

^^ , action of the final 7-nucleon system is fully taken into account via the Lorentz Integral Transform 

^D ' (LIT) method. For the calculation of the LIT we use expansions in hyperspherical harmonics (HH) 

in conjunction with the HH effective interaction (EIHH) approach. The convergence of the LIT 

y , expansion is discussed in detail. The calculated cross section agrees quite well with the available 



c 



X 



experimental data, which cover an energy range from threshold up to 100 MeV. 
PACS numbers: 21.45. +v, 24.30.Cz, 25.20.Dc, 31.15.Ja 



The electromagnetic response of an A-hodj nucleus is a basic property in nuclear physics. 
It contains important information about the dynamical structure of the system. For this 
reason microscopic calculations are needed in order to investigate the details of the reaction 
mechanism and the underlying dynamics. Traditionally such studies could only be made for 
few-body systems with A < 3. The explicit calculation of final state wave functions consti- 
tutes the main limitation of standard approaches. On the other hand these difficulties can be 

n 

avoided in the Lorentz Integral Transform (LIT) method [1] , where only a bound-state-like 
problem has to be solved although the full final state interaction (FSI) is rigorously taken 
into account. First applications of the method were made for electromagnetic responses of 
"^He 211 . Enormous progress has also been achieved in the calculation of bound-state wave 

ri 

functions of systems with A > 4 (No Core Shell Model [3|] and Effective Interaction Hyper- 
spherical Harmonics (EIHH) approach 4]). The novel methods of LIT and EIHH combined 



)led us to carry out calculations of elec- 
. Therefore, systems with 4 < A < 12 



with the help of modern computational resources ena 

tromagnetic cross sections of nuclei with A = 6 p, IS|' 

are now considered an interesting playground for establishing a bridge between few- and 

many-body physics and especially for testing many-body approximations. 

In this paper we present the first microscopic calculation of a photoabsorption reaction on 
a nucleus with A=7, namely ^Li. In this calculation we use as nucleon-nucleon interaction 
the Argonne potential AV4' [7[, which is a semirealistic central interaction that contains S- 
and P-wave forces. In recent work [6] we found that the P-wave interaction has a significant 
infiuence on the photodisintegration cross section of P-shell nuclei, like ^He and ^Li, and leads 
to a considerably better agreement with experimental data than a central S-wave interaction 
only. Thus we have chosen this potential model for the present ab initio calculation of the 
^Li total photoabsorption cross section. The calculational procedure is described in detail 
in p, la]. In the following we will briefiy review the main steps only. 

The inclusive unpolarized response function R^u) is defined by the total photoabsorption 
cross section 

a{uj) = 4:7c^aujR{uj) , (1) 

where a denotes the fine structure constant and cu the photon energy. In the LIT method [l| 
one obtains R{uj) from the inversion of the integral transform 

C{an, aj) = f du- ^^^] ^ , . (2) 



In the unretarded dipole approximation one has 



X \m{{m,{M^,M^)\b\- i -D,|vI/o(Mo^Mo^))}, (3) 

H — CT/j — lOi 

where D^ denotes the dipole operator, H the nuclear Hamiltonian, and Lq and 5*0 are the 
angular momentum and spin of the ground state |\l/o(M(^, M^)), with projections M^ and 
M^ ^ respectively (since we work with central forces, Lq and So are good quantum numbers). 
The transform C-iaR, oi) is evaluated by inserting a complete set of projection operators 
Y.cM-M^ |vl>c(M^M^))(vl/e(M^M^)| = Y^cm^m^Pcm^M^^ ^here C = {L,S,^,T^7r} 
stands for the quantum numbers characterizing the channels (angular momentum, spin, 
isospin and its projection, and parity, respectively), M^ and M^ are third components of 
angular momentum and spin. In the sum only the channels allowed by the dipole selection 
rules need to be considered and since the dipole operator does not depend on spin, we do not 
need to average over the initial spin projections (M,f ), neither to sum over M^ . Therefore 
one obtains 



ZLo + 1 a/ ^ H - an-tai 



^i^R^^i) = 7rr-T-2^\{'^cm'^o)rim{{^c\-E—: —\^c)}, 

c 



I*,) . , ^^^-l^°> , (4) 



where 



To simplify the notation, in the last two equations we have also omitted the dependences 
on M^ and M^, as will be done in the following. 

n 

In order to use the Lanczos algorithm (see Ref. [Hi) it is convenient to write C{aR, o"/) = 
YliC-^c^Ci where Afc is connected to the norm of Pc-Dzl^o); 
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^-— Y. |(vi/^|D,|vl/o)p = ^^i— Y. {'^o\D.PcD.\^o) (5) 
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ILn + ± — ^-i^n , - 



and 



Fc = -Im{(*clT7 ^ — l^c)} (6) 
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is evaluated as continuous fraction via the Lanczos coefficients. This means that the LIT of 
the total response is the sum of the Lorentz transforms of the individual channels 



C{aR,ai) = '^Cc{crR,cri) 



(7) 



c 



where for every channel one has Cc{o'r, cti) = -^c-^c- 

In the case of ^Li one has total angular momentum and parity J'^ = | and isospin 
To = I with projection Tq = — |. Using central forces, one has ground state orbital angular 
momentum Lq = 1 and spin Sq = ^. There are six different channels allowed by the dipole 
selection rules corresponding to angular momentum L = Lq — 1,Lq,Lq + 1, spin S = Sq and 
isospin T = To,To + 1 with isospin projection conserved T^ = Tq. In Table I we show the 
good quantum numbers of these channels. 

For every channel the LIT is calculated by expanding |\E'o) and |^l/c) in terms of the 7- 
body anti-symmetrized hyperspherical harmonics (HH) up to a maximum HH hyperangular 
momentum y, i.e. K^^^ for the ground state |\E'o) and K^ax for the channels I^E'c)- The 
convergence of the HH expansion is improved by using the EIHH approach. Our procedure 
consists in fixing first K'^^^ such that convergence for the binding energy is reached, and 
then we study the behavior of the LIT with increasing Kmax- As already pointed out in p, 
the rate of convergence for a potential that includes a central P-wave interaction, like the 
AV4', can be slower than that of a purely central S-wave force. Nevertheless, in the case of 
^Li we have obtained a satisfactory convergence in terms of hyperangular momentum, both 
for ground state energy and LIT. For the ground state of ^Li good convergence is reached 
with K^^^ = 9 with a binding energy of 45.28 MeV. Further increase to K'^^^ = 11 leads 
only to a small change of the binding energy by 0.05 MeV. Because of the dipole selection 
rule K = Kq ± 1 between states with hyperangular momenta K and Kq, an expansion of 



TABLE II: Reduction of basis states for channel C = 3 for K^ax = 10. Each symmetry is 
represented by a Young tableau. The total number of states for each symmetry is listed as iV*^™, 
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hsts the actual number of states in the calculation. 
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the ground state up to a certain K^^^ implies that in \^c) only states with hyperangular 
momentum K < K^ax = -^mas + ^ contribute to the LIT. Thus it is expected that for 
sufficiently high K^^^ a further increase of K^ax beyond -K'maa; + 1 "^ill '^ot result in a 
significant change. In this case a check of the convergence with respect to -ft'maa; only will 
be sufficient. 

In the present work, the best calculation of the LIT corresponds to a final state with 
Kmax = 10. In this case the number of HH-basis states to be included becomes quite large, 
especially for the channels with Lq = 1 and 2. For example, already for channel C = 3 
and Kmax = 10 one has 6348 hyperspherical states. This number has to be multiplied by 
the number of hyperradial states, about 30, to obtain the total number of states needed 
in the expansion. Therefore, it is desirable to discard those HH states which give only 
negligible contributions to the LIT. To this end we have studied the importance of the HH 
states according to their spatial symmetry and found that quite a few of them can be safely 
neglected (see Table II). We have checked this approximation by performing calculations 
with the complete set of states for lower values of Kmax (6 and 8) and compared the results 
with those using a truncated set. Whenever the differences between results with the reduced 
and the full basis were negligible (below 0.5%) we concluded that the omitted states were 
negligible also for higher Kmax- In this way we accomplished for Kmax = 10 a sizable 
reduction from A^ = 190 440 to A^ = 111 900 basis functions. In an analogous way we carried 
out the calculations for all the other channels. The estimated error introduced by these 
truncations is of the order of 0.5%. 

We are also able to check the error introduced by the symmetry truncation in a second 
way. In fact a good check of the quality of the calculation is obtained by considering the 



sum over the norms Mc defined in (jS)). Using completeness one finds 

1 



T,^'-jihri ^ i<*°i^'i*°> 



C,M'^,M^ 



2Lo + 1 ^ 



i^olDlDzl'^o) 



(8) 



where the last expression is nothing else than the mean expectation value of the operator 
D\Dz in the ground state, that can be easily calculated (see Ref. p]). With respect to 
Eq. (jHI) we obtained 1.877 [fm^] for the ground state expectation value with i^^^^. = 9, while 
using a symmetry truncated expansion for \'^c) up to K^ax = 10 we get 1.871 [fm^]. The 
small difference of 0.3% reflects the small error introduced by the symmetry truncation. 
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FIG. 1: Upper panel (a): Relative change iZ^"^''"^^'' for the LIT of channel C=6 (u/ = 10 MeV). 
Lower panel (b): i?^'""^ for the same channel (see text for details). 



As next point we address the quality of convergence of the LIT with respect to the HH 



expansion. For this reason, we first introduce the notation L 



K° (±) 



c 



for the LIT calculated 



with an expansion up to -ft'^Qj. for the ground state and up to K^^^ = K^^^ ± 1 for \^c) 



with a corresponding increase of K^ax while R^j"" represents the relative change for 



respectively. That means C^f^ represents the LIT of channel C calculated with a max- 
imal hyperangular momentum for those states which can be reached by a dipole transition 
from the ground state, whereas for Lq^ the expansion of channel C is taken up to one 

step below this value. Then we define 

Rc-'^' = ^ -j^, X 100, (9) 

'-c 

where Rq"" represents the relative percentage change for increasing K^^^^ by 2 together 

with a corresponding increase of K^ax while 

increasing K'^^^ by 2 but keeping K^ax fixed. 

In Fig. 1 (a) we show the convergence pattern of R(j"" for one of the more important 
of the six C channels, namely C = 6. One readily sees that increasing K^^^ results in a 
considerable reduction of R^ "^ . For example, changing K^^^ from 7 to 9 yields a i?g 
between 5 and 15 %, where the lower value refers to the region of the maximum of the 
LIT (compare Fig. 2), which however is still not negligible. It is very likely, that a further 
increase to K^^^ = 11 would lead to a significant improvement of the convergence. Although 
the calculation of the ground state for K^^^ = 11 can be performed, the corresponding 
calculation of the final state with K^ax = 12 is beyond our present technical capabilities 
whereas it can be done for K^ax = 10, i.e. C^ can be calculated. Thus the question now 
is, how much is the convergence improved in going from C^ to C(j . To this end we 
show R(j"'°-''^~' in Fig. 1 (b) for K^^^ = 7, 9, 11. Again one notes a considerable decrease of 
j^^maxi-) ^[j^\^ increasing K^^^ and for the highest value K^^^ = 11 of the order of 5 % in 
the region a^ < 60 MeV where the LIT is more sizable. Furthermore, comparing panels (a) 
and (b) of Fig. 1 one can see that for a given K'^^^ the relative change Rq "^ is smaller 
than that of Rq """"" . Therefore, we expect that a further increase to £g would lead to 
a relative change with respect to £g smaller than 5%. Since the other channels behave 
similarly, we estimate a total uncertainty of about 5%. 

Now we turn to the results for photodisintegration. In order to obtain the total pho- 
toabsorption cross section cr{u!), we first invert the LITs of the response functions for the 



various channels (for details see ll|). Then one obtains the contribution to the cross sec- 
tion of each channel using the relation (^. Summing up all contributions finally yields the 
total cross section a{uj). In Fig. 2 we show both the calculated total cross section and the 
contributions of the separate channels. The energetically lowest open channel is the T = \ 
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FIG. 2: Contribution of various channels to the total cross section. Panels (a) and (b) show 
the separate contributions of the different channels and their sum for T = 1/2 and T = 3/2, 
respectively. Panel (c) shows again the T = 1/2 and T = 3/2 contributions and the total cross 
section. 

associated with the reaction ^Li +7 — > "^He +t with a threshold of 2.47 MeV (the theoretical 
threshold obtained with the AV4' potential is 4.70 MeV), whereas the lowest open T = | 
channel corresponds to ^Li +7 -^ ^He + p, whose threshold is 9.975 MeV (the theoretical 
value is 12.41 MeV). In Fig. 2 the theoretical thresholds have been used. 

For the T = 1/2-channels in Fig. 2 (a) one readily sees that by far the largest contribution 
comes from channel 3 rising steeply above threshold, reaching a maximum around 17 MeV 
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FIG. 3: Comparison of the theoretical photoabsorption cross section calculated with AV4' potential 
with experimental data from [12:] . 

and falling off only slowly. Channel 5 is the next in importance rising only slowly above 
10 MeV with a maximum near 33 MeV and becoming then comparable in size to channel 
3. Only in the very near threshold region channel 1 is dominant but then becomes much 
smaller than the other two channels. In Fig. 2 (b) the T = 3/2 channels have two dominant 
contributions, almost similar in size. Channel 4 is slightly larger showing also a steep rise at 
threshold and slow fall-off with a maximum near 23 MeV whereas the second in importance, 
channel 6, shows a slow rise and a peak around 40 MeV. Compared to these two channels, 
the remaining channel 2 appears quite marginal. In view of the two maxima of almost 
equal height with a separation by about 17 MeV the total T = 3/2-contribution exhibits a 
broader distribution than the T = 1/2-contribution with a shoulder on the low-energy side. 
The maxima of both contributions have about the same size but are separated by about 
20 MeV. Thus, the resulting total cross section in Fig. 2 (c) shows also a broad distribution 
with a steep rise right above threshold, a slight shoulder above the maximum and a slow 
fall-off at higher energies. 

This characteristic behavior is indeed exhibited by the experimental data on ^Li in Fig. 3 



IJ. Note 



where we show a comparison of the theoretical result to experimental data from 
that the theoretical cross section is shifted here from the theoretical threshold to the ex- 
perimental one. One readily notes that the gross properties of the data, steep rise, broad 
maximum and slow fall off, are very well reproduced quantitatively over the whole energy 



region by the theory. It is worthwhile to emphasize that this result is based on an ab initio 
calculation in which the complicated final state interaction of the 7N-system is rigorously 
taken into account by application of the LIT method. No adjustable parameters were used, 
the sole ingredient being the AV4' NN potential model. It remains to be seen whether the 
slight variation of the data near and above the maximum will also be found in an experi- 
ment with improved accuracy. Therefore, a new measurement of the total cross section with 
a higher precision would be very desirable. In particular, this could clarify the question 
whether a simple semi-realistic potential like the AV4' model is sufficient for an accurate 
theoretical description of this reaction or whether a more realistic nuclear force including a 
3N-force is needed. 
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